Scattering of circularly polarized light by a rotating black hole 
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We study scattering of polarized light by a rotating (Kerr) black hole of the mass M and the 
angular momentum J. In order to keep trace of the polarization dependence of photon trajectories 
one can use the following dimensionless parameter: e — ±{ijjM)~^ , where oj is the photon frequency 
and the sign + (— ) corresponds to the right (left) circular polarization. We assume that |e| <C 1 and 
use the modified geometric optics approximation developed in that is we include the first order 
in e polarization dependent terms into the eikonal equation. These corrections modify late time 
behavior of photons. We demonstrate that the photon moves along a null curve, which in the limit 
e = becomes a null geodesic. We focus on the scattering problem for polarized light. Namely, we 
consider the following problems: (i) How does the photon bending angle depend on its polarization; 
(ii) How does position of the image of a point-like source depend on its polarization; (iii) How does 
the arrival time of photons depend on their polarization. We perform the numerical calculations that 
illustrate these effects for an extremely rotating black hole and discuss their possible applications. 



PACS numbers; 04.70.Bw, 04.25.-g, 42.15.Dp 



I. INTRODUCTION 

There exists a well-known analogy between gravity 
and electromagnetism. In particular, the linearized Ein- 
stein equations can be written in the form similar to the 
Maxwell equations (see, e.g., the review [2,] and references 
therein) . Analyzing propagation of a circularly polarized 
light in a rotating frame in a flat spacetime, Mashhoon 
[3| demonstrated that there exists helicity-rotation cou- 
pling with the energy proportional to ±h{k,^), where 
k is the wave vector of a photon, and $7 is the angular 
velocity of the rotation. This effect is similar to energy 
splitting of a particle with magnetic moment in a mag- 
netic field. When such a particle moves in an inhomoge- 
neous magnetic field its trajectory depends on the spin 
orientation. One can expect that there exists a gravi- 
tational analogue of this Stern- Gerlach effect. Namely, 
in the gravitational field of a rotating body trajectories 
of circularly polarized photons depend on their polariza- 
tion. Using this analogy, Mashhoon demonstrated 
that photons of the opposite circular polarization emitted 
by a distant source, after scattering, deflect to directions 
with the separation angle. 



separation angle ■ 



\GJ 



(1) 



Here A is the wavelength of the photon, J is the angular 
momentum of the rotating body, and D is the characteris- 
tic distance from the photon to the body at the moment 
of their minimal separation. Mashhoon also arrived to 
the conclusion that in order to describe these effects in 
the geometric optics approximation, the corresponding 
corrections should be included in the eikonal approxi- 
mation. It should be emphasized that these conclusions 
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are based on the gravitational Larmor theorem Q and 
gravito-electromagnetic analogy, valid in the weak fleld 
approximation. 

One can expect that the dependence of the photon 
trajectory on its helicity might be enhanced in the case 
when the photon moves in a strong gravitational field, for 
example close to a rapidly rotating black hole. In order to 
derive this dependence, one needs to modify the standard 
geometric optics (GO) approximation and adapt it to the 
propagation of circularly polarized beams of light. The 
standard GO scheme is well known (see, e.g., [6]). In 
this approach one starts with the following anzats for 
the vector potential of the electromagnetic field in the 
Lorentz gauge: 



)• 



(2) 



One assumes that e is a small parameter, so that the 
wavelength of light A is much smaller than the other 
scales L, e = A/L ^ 1. After a substitution of this 
anzats into the Maxwell equations 



(3) 



one obtains in the lowest order approximation the follow- 
ing eikonal equation: 



(4) 



This first order partial differential equation can be solved 
by the method of characteristics by identifying V^S* — 
kfj_, where k^^ is a vector tangent to null geodesies de- 
scribing the photon trajectory. The next order in e equa- 
tions show that the vector of linear polarization is par- 
allel transported along the worldline of the photon and 
the following conservation law is satisfied: 



V^(a^F 



0. 



(5) 



Higher order corrections give expansion of the amplitude 
flp in terms of the small parameter e, while the eikonal 
function S is kept unchanged. 
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For a real scalar field the anzats similar to ([2]) is well 
defined and unique if one chose the prefactor to be real. 
For a many component field there is an ambiguity: one 
can add to S higher order 'correction' eip and, at the 
same time, change the prefactor — ?> exp{—iif)a^ (for 
discussion of the WKB and GO approximations for multi- 
component field see, e.g., Q)- A similar problem arises, 
for example, when one considers WKB approximation for 
the Dirac equation describing an electron moving in the 
external electromagnetic field (see, e.g., The re- 

cent nice paper [l^ clearly illustrates problems that arise 
when one tries to adapt WKB approximation for the de- 
scription of the Stern- Gerlach effect. It is shown that the 
correct answer is obtained when one first diagonalizes 
the corresponding Pauli equation, and then includes the 
spin corrections, which are of the first order in h, into 
the eikonal equation. Namely, these first order terms 
modify the propagation of the electron at late time (at 
far distance) and allow one to obtain the correct result. 
This result directly follows from the stationary phase ap- 
proximation for the path integral for the electron in the 
magnetic field. As a result of this approximation, one 
restores the theory of the electron as a particle with 'a 
classical spin', and the corresponding Dirac (Pauli) equa- 
tion is nothing but the first quantized theory. 

Unfortunately, a similar approach does not work for 
an electromagnetic field. The problems are basically con- 
nected with gauge invariance. A massless photon can be 
easily defined in the momentum space, while it is a non- 
trivial problem to define a localized classical photon in 
the configuration space (in the spacetime). The compre- 
hensive discussion of this problem can be found in [17| . 
see also the references therein. 

In our paper [ij , in order to study propagation of circu- 
larly polarized monochromatic electromagnetic waves in 
a stationary spacetime, we used another approach. Using 
the projector along the direction of the Killing vector, the 
Maxwell equations are reduced to (3 -I- 1)— dimensional 
form. The three-dimensional complex vectors — 
EztiH, that are linear combinations of the electric E and 
magnetic H fields, are used to describe fields with fixed 
circular polarization. In a flat spacetime are known 
as the Riemann-Silberstein vectors (see, e.g., [l^l)- It is 
easy to show that in a stationary spacetime their positive 
frequency Fourier transform obeys a single master 
equation, which describes a state with a given circular 
polarization. 

In P, we used the GO anzats for .F^, 

JP± = /±e'^*/^ . (6) 

For monochromatic waves of high frequency e ~ uj~^ is 
a small parameter. We demonstrated that there exists a 
natural prescription for the choice of the eikonal function. 
According to this prescription, the lowest order helicity 
dependent corrections should be included in , while 
the transport equation for the polarization vector (at 
least in the lowest order in e) does not depend on the 
state of polarization. In such a modified GO approach the 



eikonal equation, which determines the effective Hamilto- 
nian for light rays, includes terms proportional to e. The 
sign of this terms depends on the helicity. As a result, 
the modified Hamiltonian describes rays with worldlincs 
depending on the polarization of the photon. For short 
time asymptotics such a modification may be considered 
as a simple re-arrangement of terms in the asymptotic ex- 
pansion in the parameter e, while in the late time regime 
it results in the spatial separation of the light rays for 
photons with opposite circular polarization. For initially 
linearly polarized light at late time one obtains a mixture 
of two opposite circularly polarized photons, which are 
moving along spatially different trajectories. 



In a stationary spacetime this effect exists only when 
the gravitating body is rotating. If t is the Killing time, 
then for a rotating mass the gravitational field contains 
components gu {i = 1,2,3) which cannot be banished 
by the spatial coordinate transformation x*(a;*). 
The remarkable fact is that additional corrections to the 
3D equations of motion of the photon connected with 
its polarization can be described as the additional force 
±e curl curl g in these equations 



The goal of this paper is to study the scattering of a 
circularly polarized light by a rotating (Kerr) black hole. 
The paper is organized as follows. In Sec. II we in- 
troduce ultrastationary metric for stationary spacetime 
and discuss (3 -t- 1)— form of the polarized ray equations. 
We demonstrate that circularly polarized beams of light 
propagate along null rays, which become null geodesies 
in the limit e = (that is, for unpolarized light). We also 
describe some general results concerning time of arrival 
of polarized photons emitted at some point to a distant 
observer. In Sec. Ill we apply this formalism to the 
Kerr spacetime. In addition, we introduce the advanced 
and retarded star time coordinates which increase mono- 
tonically along the null rays and which are convenient 
for their parametrization. We present the equations for 
the propagation of circularly polarized photons in dimen- 
sionless form which is suitable for numerical integration. 
Section IV contains formulation of the scattering problem 
in the Kerr spacetime adapted to the circularly polarized 
photons. In Sec. V we study how polarization affects 
bending angles of scattered photons, and in Sec. VI we 
study how position of the image of a photon and its time 
of arrival depend on polarization. We discuss the results 
in Sec. VII. Details of the numerical computations are 
given in the Appendix. 



In the present paper we use the system of units where 
G — c — 1 and the sign conventions for the metric and 
other geometrical quantities adopted in the book Q. 
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II. CIRCULARLY POLARIZED LIGHT 
PROPAGATION IN A STATIONARY 
SPACETIME 

A. Effective equation for polarized photons 

Let us consider the following equation in a curved 
spacetime with the metric (7^^,: 



P ^eg^PF,^—, (7) 



d\ 



where D/dX is the covariant derivative, P is a scalar func- 
tion, e is a small parameter, and Fp^ — f [pi/j • We shall 
demonstrate later that the equation for propagation of 
circularly polarized light in a stationary spacetime is of 
this form. Let us at first discuss general properties of this 
equation. Let us multiply Eq. JT]) by g^^^.dx'^ /dX, then 
one has 



(8) 



dV dx^' dx" 

This relation implies that if for some value of A the func- 
tion V vanishes, then V = along the complete world- 
line, and hence this line is null. The reparametrization 



dX = ^dX , 



(9) 



where /3 is a scalar function, preserves the form of Eq. 
([7|), and one has 

P^+e(3g^''Fp,^, (10) 



dX^ 



dX 



dX 



dln/3 
dX 



(11) 



Thus, by proper rescaling one can always put P = 0. 
We call such a parametrization, defined up to a linear 
rescaling, an affine parametrization. For e = 0, Eq. ([7]) 
is a geodesic equation. For null geodesies the canonical 
parameter A is a usual affine parameter. 

Consider now a conformal transformation of the met- 



B. Stationary spacetime and ultrastationary metric 

To study how polarization affects the light propagation 
in a stationary spacetime, it is convenient to use what 
we call ultrastationary form of the metric. Let us discuss 
the related formalism. We present the metric ds^ of a 
stationary spacetime in the form 

ds^ = gf.udx'^dx" = hdS"^ , h = -gu = -(.ft) 1 (15) 
dS^ = -{dt-g,dx'f + dl\ g^^gu/h, (16) 
dp = 7y dx^ dx^ , 7ij- = h^^gij + g^gj . (17) 

Here ^^^^ = (5^ is a timelike at infinity Killing vector. For 

gi = the metric dS^ is ultrastatic. We call the metric 
dS^ ultrastationary. 

Source-free Maxwell equations in 4D are conformally 
invariant. For this reason, in order to study their solu- 
tions, one can use either the original spacetime metric 
ds^ or its ultrastationary version ^5*^. Similarly, a null 
ray in the physical metric ds'^ is at the same time a null 
ray in the conformal metric dS^. Moreover, after the 
conformal transformation, a null geodesic remains a null 
geodesic. However, its affine parametrization depends on 
the conformal factor. In our discussion we shall use the 
ultrastationary form of the metric, that simplifies essen- 
tially many relations. 



C. (3 + 1) form of the equations 

Because of the conformal invariance of a null ray so- 
lutions to the effective equation (fT4)) . it is sufficient to 
study their properties in the ultrastationary spacetime. 
Consider a curve x^{X) and denote by = x^ = dx^^ /dX 
its tangent vector. We also use the notation 



(...)• =u^V^(...) 



(18) 



where is the covariant derivative in the ultrastation- 
ary metric dS^ . In these notations the effective equation 



gt,u = ^'^gtj.v ■ (12) 
For null worldlines Eq. ([T]) is transformed as follows: 



rflnOA dxf" , „„ dx" , , 



By comparing this result with Eq. ((T0| it is easy to see 
that the canonical form of Eq. ((T]), 



D'^xi' 
dX^ 



dr" 
eg tp, , 



(14) 



is conformal invariant, provided under this transforma- 
tion F^i, = P^^, and the conformal transformation is 
accompanied by the reparametrization ^ with /3 = Jl^ 



^ It is interesting to notice that Eq. I I14II can be obtained as a 



limiting case of the equation of motion of a charged particle of 
the mass m and the electric charge e. This equation is of the 
form 

"^^ = '=^'^^^' 
where D/dr is the covariant derivative corresponding to the met- 
ric gfiv. Denote by A a parameter connected with the proper time 
T as follows: 

T = mX . 

Then, this equation takes the form 

where the over dot stands for the covariant derivative D/d\. It 
is easy to check that 

= = —m? . 

The above equations allow for the ultrarelativistic limit m — > 0, 
and in this limit one obtains Eq. II14I I with e = e. 
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(iMl) takes the form 



(19) 



Let us write this equation in a (3 + 1) form constructed 
by projecting it along the Kilhng vector ^^j. We notice 
that the only non- vanishing components of the Christoffel 
symbols with the only one index t are 



Here the covariant derivative D/dl is defined with respect 
to the metric jij. In what follows, all the operations 
involving the Latin indices j, . . . are performed by using 
the metric jij and its inverse 7*^ . In particular, the vector 
and the scalar products of vectors a and b are given by 



(a, b) = a'-jtjV , [axbY = e^^'^ajhk , 



ijk 



r 



tij 



(20) 



Denote by Qijk the ChristofFell symbols for the 3D metric 
7ij. Then one has 



+ 9k{9js - 9ij) - 9i{9],k + 9k,j)] 



(21) 



Denote = g^^x'^ , where 5^1/ is the ultrastationary met- 
ric. Then one has 



Denote 



'9r 



d^x' 
d)^ 



dx^ dx^ 



-IT ■ (22) 



respectively. Here e*^*^ is the three-dimensional com- 
pletely antisymmetric Levi-Civita symbol normalized by 
the condition e^^^ = 1. 

It was shown in that the equations for propaga- 
tion of circularly polarized monochromatic light of the 
frequency w have the form ([^^ -([50 ]) with 



F,, - 2<i>[,,,] , $^ = -(curlg)^ - [V x gY = e'^'^g.^k ■ 



F ~ curl curl g . 



(32) 



One also has e = ±(2a;M)~^, and the signs (— ) and 
(-I-) define the left and the right circular polarization, 
respectively. In what follows, we use the subscript (0) for 
a null geodesic, while the subscripts (— ) and (-I-) stand 
for the left and right polarization of photons, respectively. 



dX 



dx' 



(23) 



D. Time delay and Fermat principle 



then Eq. 



can be written in the form 



dU 
"dX 



Using this relation and Eq. (|23| one obtains 



lij- 



dA2 



, .dx^ dU 



(24) 



(25) 



where D/dX is the covariant derivative corresponding to 
the 3D metric 7ij. 

Let us consider a special case of the field when 
Fu = 0. Then the 4D equation ((M]) has the following 3D 
form: 



dU 
'dX 
dt 



0, lij- 



D^x^ 



dA2 



dx^ 
'~dX 



dx^ 

^ in 



'U{gi,j - gj,i) + eF,. 



(26) 

(27) 
(28) 



Here U is the constant of motion. By the rescaling A — > 
AX (A — const) one can always put U ~ 1. In what 
follows, we shall use this choice. 

Denote x = {x^,x'^,x^), then Eqs. (I26l)-(l28l) can be 
written in the following 3D form: 



~dP 



dx 

— X (curlg + eF) 



dt _ 



dx 



(29) 
(30) 



Suppose r is a timelike worldline of an observer and 
t is the proper time along it. Suppose at some moment 
of time to the observer registered light emitted at some 
point e. Consider also null rays emitted from the same 
point e and obeying Eq. (HH). We assume that for each 
of the value of the parameter e from a small interval 
there exists a null curve intersecting F. Let us denote 
by t{e) the time of its arrival to F. Thus, one has a 
one-parametric family of null rays emitted at e and in- 
tersecting F. It was proved in the papers p^l - l2l| that 
{dt{e)/de)\e=o = 0, and if there is no conjugated points 
along the null geodesic, then t{0) is minimum of t{e). In 
other words, in the general case the null geodesic arrives 
to the observer F earlier than other null curves. 

This general result concerning the time delay of null 
rays takes a simple form in a stationary spacetime. 
Namely, let Fg and Fq be two timelike Killing trajec- 
tories, and Xe and Xq be two spatial points representing 
these trajectories in the 3D space. Let a photon emitted 
at Xe at the time te is registered by an observer at Xo at 
the time to- The relation (jSHl) shows that 



t{e) = te+ {dl + gidx') 



(33) 



It is possible to show (see, e.g., [2^l23|) that for a one- 
parametric family of null curves close to a null geodesic, 
the latter gives the extremum of the arrival time t{e). 
In the general case, this extremality means that a pho- 
ton moving along such a geodesic arrives to the point of 
observation earlier than any other photon which follow 
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non-geodesic motion. Possible exception is the case when 
there exist more than one null geodesies connecting the 
emitter and the observer. This is nothing but the Fermat 
principle for light rays in a stationary gravitational field. 



III. PROPAGATION OF CIRCULARLY 
POLARIZED LIGHT IN THE KERR SPACETIME 



For photons the functions R and 8 are 



R=[{r 



'')E - aL]^ - A[Q + (L- aE'f] , 



e = Q + a^E^ cos^ 0^ cot^ 



(46) 



Important property of a null geodesies of a Kerr space- 
time is that it has no more than one radial turning point 
( see, e.g., gl). 



A. The Kerr spacetime 

The Kerr metric in the Boyer-Lindquist coordinates 



x" = (t, r, 9, (/)) can be presented in the form 
ds^ ^ h[-~{dt~ gidx'f +df] , 

where 



h = l- 



2Mr 



dp = "fijdx^dx-' 
and 



9i = - 

~ ~Kh 



2aMr 



■ sm 



dr^ 



-dO^ 



A sin^ e 



S = + cos^ ( 



A = - 2Mr 



(34) 

(35) 
(36) 

(37) 



This metric describes the gravitational field of a black 
hole of the mass M and the angular momentum J = aAI 
(a > 0, a/M < 1). 



The Kerr metric has two Killing vectors ^^^^ and ^; 



and the second rank Killing tensor i^^^ iK(p,^.a) = 0)- 
For a geodesic line x^(A) with an affine parameter A 
(m'^ = dx^/dX) the corresponding conserved quantities 
are (see, e.g., ^) 



E : 



, I dt dcj) 



dt 



h 9r 



K0) " \'"^d\ 
Q = K-{L- aEf , K ^ K°'^Uo,up 



dX 



(38) 

,(39) 
(40) 



Here E is photon energy, L is photon azimuthal angular 
momentum, and Q is the Carter constant. These inte- 
grals of motion make the geodesic equation 



D'^xf' 







(41) 



completely integrable. In particular, this means that the 
geodesic equation can be written in the following first 
order form: 



dt 

dX 

dr 

dX 
dd 

dX 



(.^+a^)^^-2aM.L_^,^.^,^^ ^^^^ 



E — 

dX 



aj-Vli , <Jr — ±1 



creV6 , fTfj = ±1 , 
L 2aMrE - a'^L 



(43) 
(44) 
(45) 



B. Advanced and retarded star time 

In the study of the Kerr spacetime one often uses the 
retarded and advanced time coordinates 



u = t — 



V ^ t + r^, 



(r^ + a'^)dr 
A 



(47) 



They are well-known generalization of the outgoing 
and ingoing Eddington-Finkelstein coordinates of the 
Schwarzschild spacetime (see, e.g., Q). The latter are 
null coordinates, that is, surfaces u = const and v = 
const are null. However, as it was demonstrated in f25|. 
the coordinates u and v in the Kerr spacetime do not 
have this property. Namely, in the presence of rotation, 
surfaces u — const and v — const are timelike, and they 
become null only at infinity. For this reason, one can- 
not guarantee that u and v are monotonically increasing 
functions along any null curve. 

However, it is possible to introduce new coordinates 
which are almost everywhere timelike and which can be 
used as monotinic parameters along null curves. We de- 
fine these coordinates as follows: 



W± 



H 



{r* + 2a^Mr 



(48) 
(49) 



It is easy to see that 



V — = W+ — u = r,t~ri,, dr^ — dr^, = — dr . 

(50) 

To distinguish from the tortoise coordinate , we call 
it radial star coordinate. 

Expanding the last expression in Eq. (|50p at r — > oo 
one obtains 

/ 3a2 la^Aj 
dr* - dr* — 1 - t ^ - — ^ 



+ ... dr. (51) 



This means that at a large distance the difference be- 
tween W_ and and the corresponding values of v 
and u is of the order of a^/r, so that the new coordinates 
W- and W+ are also asymptotically null. 

Let us show that a surface W± — const is always non- 
timelike, so that W± can be used as a time coordinate. 
This property differs W± from u and v. Indeed, the gra- 
dient covector W±;^ in the Boyer-Lindquist coordinates 
reads 



1,T^,0,0 



(52) 
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Using the expressions for the components of g'^'^, 



then the Kerr metric (l34l) takes the form 



2^2 



„ a^Asin^ 6* - (r^ + a^) 
9 = 

one gets 



AS 



A 



> 0. 



(53) 



(54) 



For a rotating black hole the equality is valid only when 
6 — tt/2. In other words, the gradient V1V± is always 
non-spacclikc. It is timelike everywhere outside the equa- 
torial plane and is null on this plane. We call the coor- 
dinates W- and W+ constructed by using the radial star 
coordinate r^, advanced and retarded star time, respec- 
tively. 

Suppose ^''(A) is a null curve in the Kerr spacetime, 
then the change of the coordinate W± along this curve is 



dW± dW± dx" 



dX 



dx°' d\ 



> 0. 



(55) 



The equality is valid only if a null curve crosses the equa- 
torial plane, where W±;^ is null, and at that point its 
tangent vector dx'^/dX is colinear to 

W^\g=^/2^-{r^^)-HH^,±HA,0,2aMr). (56) 

If one excludes these special null curves, one can use both 
the star time coordinates W+ and W- as monotonically 
increasing parameters along null curves. The Kerr metric 
(p4|) in the coordinates x" = {W±,r,6,(j)) remains the 
same with the only change 



dt — gidx^ — ^ dW± — g^dx^ 



2aMr 



■ sm 



(57) 
(58) 



C. Dimensionless form of the equations 

The mass parameter M in the Kerr metric determines 
the scale. In what follows, it is convenient to use the di- 
mensionless form of the dynamical equations given above. 
Instead of the coordinates {t, r, 9, (f) and W± we shall use 
the coordinates (t, cc, y, <j)) and w± defined as follows: 



t 



M 



- ^± + /■ ^ w 



(59) 



where 

rj={l + 2a^x^ + a^x'^Y/'^, 5 = I - 2x + a^x'^ . (60) 
Let us denote 

a = ^ , di^ = M-^dl^ , (61) 



Here 



ds^ = M^hdS^ , 

dS^ = -{dr-jidz^f + df. 

2ax{l - 2/2) 



(62) 
(63) 



2x 
a 



It 



(g - 2x) 



2 2 2 

a X y . 



h = 1 

cr = 1 

are spatial coordinates, = (a;, y,0), and 
where 



Ixxdx^ + lyydy"^ + j^^dcj}'^ , 



x^5{(j - 2x) 
a;2((7-2x)2 



a;2(cr-2x)(l-?/2) 



7 = dct(7ij) = 



x^{(j~2xY 



(64) 
(65) 

(66) 
(67) 



In these coordinates the horizon x = xh and the ergo- 



sphere x = x 
the relations 



of a Kerr black hole are determined by 



1 



XH 



^erg 



1 



(68) 



The equation for propagation of the circularly polar- 
ized photons in the ultrastationary spacetime (|63p , writ- 
ten in the (3 + 1) form and expressed in the dimensional 
coordinates has the following form [cf. Eqs. (|29l) . (|30l) . 
and ((321) ]: 



dr 
dl 



dz 



1, 



dz 



fs = curl 7 + ecurl curl 7 , e = ± 



1 



2ujM 



(69) 
(70) 
(71) 



As we explained above, w is the frequency of a photon at 
infinity, and the sign (±) defines its helicity, which cor- 
responds to the right (-I-) and to the left (— ) circularly 
polarization of an electromagnetic wave. We always as- 
sume that |e| ^ 1. In the Kerr metric the no n- vanishing 
components of curl 7 and curl curl 7 are given by 



(curl 7)3. 
(curl 7)2, 
(curl curl 7)0 



Aaxy 
a{<T - 2x) ' 
2ax'^{l — a'^x'^y'^) 

a{cr — 2x) 
4q!x'*^(1 - y2) 



(72) 
(73) 
(74) 



cr(cr- 2a;)2 

For £ = 0, that is, for null geodesies of the Kerr metric. 
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Eqs. ((69|) - (|7ip can be written in the first order form 

^ = a,^{cT - 2x)Vx , a, = ±1, (75) 
at (T^ 

^ = '^?^J(^-2a;)\^, a, = ±1, (76) 



^ ^ (p - g^) 



(77) 
(78) 



Here 



X = [l~ ax^ip ~ a)f - x^S[q + [p - af] , 
3^=(g + aV)(l-y')-pV, (79) 
where p and g are dimensionless conserved quantities, 



EM ' ^ E^M^' 



Q 



(80) 



IV. 



SCATTERING PROBLEM IN THE KERR 
SPACETIME 

A. Scattering data 



Our goal is to study propagation of a circularly polar- 
ized photon in the Kerr spacetime. Its trajectory is a 
solution of Eqs. (|M)) - ([7T|) . We always assume that the 
coefficient e which controls the relative strength of the 
polarization corrections is small, and in the dynamical 
equations consider the terms with non-vanishing e as a 
perturbation. Let us focus on the scattering problem. 
We assume that a photon begins its motion at the past 
null infinity J^^ , propagates towards a Kerr black hole, 
gets scattered, and finally escapes to the future null in- 
finity j/'*' . Certainly, some of the incoming photons can 
be captured by the black hole and do not reach J'^. We 
do not consider them in this paper. Moreover, we use 
the (3-1- 1)— formalism (see, ^|), which is well defined in 
the domain where the Killing vector is timelike. Hence, 
we assume that the photons propagate outside the black 
hole ergosphere. 

The spin-optical interaction described by the effective 
force /e is proportional to curl curl g. It grows when a 
photon approaches the black hole and reaches its max- 
imum near its radial turning point. Far away from the 
Kerr black hole this extra force rapidly falls down as 



I curl curl g \ 



(81) 



Hence, the polarization dependent corrections to a pho- 
ton trajectory become negligible in the asymptotic do- 
main where r — >■ oo. In other words, photons of the left 
and right helicity practically follow null geodesies in the 
asymptotic region, and they can be uniquely specified by 
the same 'scattering data' as unpolarized photons. 



Let us remind these 'scattering data'. To be concrete, 
let us consider the future null infinity . This is a three- 
dimensional surface which can be parameterized by the 
angular coordinates and 0, which specify the direction 
of motion of an outgoing photon and the moment of its 
arrival to J7+, W+ — Mw+. It should be emphasized 
that the star time coordinates w± are defined up to a 
constant. To fix this constant, one needs to chose the 
initial radial coordinate in the integral for in (|48|) . We 
shall discuss this point in Sec. V. In fact, a set of parallel 
moving photons arrive at the same point of ^7"*". The 
'members' of this set can be parametrized by two addi- 
tional parameters, ^ and C, called the impact parameters. 

At the large distance solutions of the photon equa- 
tion of motion in the Kerr metric are of the form [cf. Eqs. 



dr 



L 



dt Er'^ sin 



(82) 
(83) 



d9 



1 



^ - ±^JQ + E^a^ cos2 6*0 - cot^ Oo . (84) 
dt Er^ ^ 

The angles of displacement of the photon in the (j) and 9 
directions are r ^ sin 9 od(j)/dt and rod9/dt, respectively. 
These angles decrease as r~^. Multiplying these quanti- 
ties by To one obtains the impact parameters (M^, M(^). 
The dimensionless impact parameters (^, C) are related 
to the asymptotic integrals of motion (p, q) as follows: 



p 



smf 



, Co = a-0\ q + a'^ cos2 



p2 cot 



(85) 

The impact vectors p — (^, C) form a two-dimensional 
impact plane shown in Fig. [T] In what follows, it is 
convenient to use the polar coordinates {p, x) on this 
plane. To summarize, a photon arriving at is 
uniquely specified by the five dimensionless parameters 
Vo = {w-o,yo = cos 6*0,00, Co, Co}- We call these pa- 
rameters the scattering data for an outgoing photon. A 
similar set = {w+e^ye = cos9e,(f>e, Pe,Xe} forms the 
scattering data for an incoming photon. 

The scattering problem for photons in the Kerr space- 
time can be formulated as follows. Determine I?e for a 
photon which begins its motion with the parameters I?o ■ 
One may consider a map 



(86) 



as a classical analogue of an S*— matrix. In the absence 
of spin-optical effects (for null geodesies) the scattering 
problem can be reduced to elliptic integrals. In the pres- 
ence of spin-optical corrections, when the equations for 
a null ray contain the parameter e, this map is slightly 
modified and becomes 'l'(e). We would like to study this 
modified spin-dependent map as a function of the pa- 
rameter e. In accordance with our main assumption we 
restrict ourselves to effects linear in e. Note that in the 
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FIG. 1. Impact plane. The direction of ^ is in the positive 
direction of (j) (direction of the black hole rotation), and the 
direction of is in the positive direction of 9. A circle of 
p = const corresponds to a constant asymptotic value of a 
photon total angular momentum. 



expansion of ^'(e) with respect to e, the zero and the 
first order in e terms are nontrivial functions of the coor- 
dinates {w±,x,y,(f>), which are defined by a solution to 
the dynamical equations (|M | - ((7T|) . 



1+ 




J- 

FIG. 2. Diagram illustrating scattering of photons by a Kerr 
black hole. The angular coordinate 6 is suppressed. A pho- 
ton trajectory is shown by the curve connecting the points e 
(emitter) at J'~ and o (observer) at J'^ . 



The scattering of photons by a Kerr black hole is 
schematically illustrated in Fig. [51 If the black hole 
is absent, i.e., if the spacetime is Minkowski one, then a 
photon trajectory from the emitter to the observer is a 
straight line. In this case, one can easily define the fi- 
nal angular coordinates ((/)o',0o') at J^~^ in terms of the 



initial ones {4>e,de) at J' as follows: 

do' =TT-de, (l>o' ^ TT + (f)e ■ (87) 

In the presence of a black hole, a photon trajectory will 
deviate from the straight line, so that (j)o — ipe > ■ To 
measure such a deviation, one can define the bending 
angles 9o — {tt — Oe) and (f>o — {tt + cj)e) > 0. In the 
absence of spin-optical interaction, photons propagate 
along null geodesies of the Kerr spacetime. In this case, 
the bending angles can be calculated analytically. In the 
case when the angle 9 has only one minimum or maxi- 
mum, the bending angles were calculated in [26j up to 
the third-order terms (M/rmin)^, {a/rminY, where rmm 
corresponds to the radial turning point of a null geodesic. 

B. Impact parameters and asymptotic integrals of 
motion 

There exist two different ways how to specify a null 
curve reaching at the given values {'w±,y, (/)): (i) ei- 
ther to give its impact parameters, or (ii) to fix its asymp- 
totic integrals of motion p and q. Let us discuss this issue 
in detail. Equation ([55)) establishes the relations between 
these two different sets of data. Near J'~ one can use the 
Cartesian coordinates {t, X, Y, Z) . In these coordinates 
an incoming photon of the energy E has the following 
4-momentum defined at J^~: 

p = E{— sin^e cos(/)e, — sin 6*2 sm(t>e^ — cosOe) , (88) 

while the spatial dimensionless vector p corresponding to 
the impact parameters (Ce,Ce) has the coordinates 

PX = -ie sin (t>f, -I- Ce cos 9e COS 0e , 
PY = COS 0e -f Ce COS Of, sin (/)e , 

pz = -CeSin6le. (89) 
Calculating the photon total angular momentum 

L^M[pxp] , (90) 

we obtain 

p = pcosxsin0e, p=^WTQ^\L\/{EM). (91) 

It is possible to show (see, e.g., [23|) that in the asymptot- 
ically flat region the Carter constant Q and the asymp- 
totic integrals of motion are related as follows: 

=Q + Ll + a^p% . (92) 

Using the relations ([91]) and ([92]) one obtains 

q = - p^ cos^ X sin^ 9^ - c? cos^ 9e ■ (93) 

Equations ([9T|) and ([93]) are equivalent to Eq. ([85]) . For 
s = and the given initial point {w_,,,9,,,(j),,) these in- 
tegrals of motion uniquely determine a null geodesic. 
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For £ 7^ 0, p is still an integral of motion, while q is 
only asymptotic invariant. However, [p, q) defined at 
(it;_e, ^e, 0e) Uniquely determine a null ray solution of 
the perturbed equations ((M1) - ([7T|) . 

Let us fix the star time W-e of a photon emission and 
the spherical coordinates {9e,(j)e) of the emission point. 
In what follows, we discuss the following three problems: 

f . How does the photon bending angle depend on its 
polarization? 

2. How does position of the image of a photon arriv- 
ing to an observer at {9o,(j)o) depend on its polar- 
ization? 

3. How does the arrival time of such photons depend 
on their polarization? 



spherical coordinates which are evidently behave badly 
near the axis of symmetry 9 — 0, tt. In particular, for 
null rays passing close to the axis of symmetry, a rela- 
tively large value of A0o does not mean that the bending 
angle is large. A better measure for the polarization shift 
of the bending angles is the quantity 

Aa = [(A0o)2 + sin2 0o(A0„)2] . (95) 

This is the invariant distance on a unit (celestial) sphere 
between its two close points, divided by |e| . To keep trace 
of the direction of the displacement we shall use the angle 
/3 between the displacement vector and the (/)— direction, 

/\9 

tan/3 = . ° , ■ (96) 
sm 9o A(po 



V. POLARIZATION SHIFT OF THE BENDING 
ANGLES 

A. Shift parameters 

Let us discuss the first problem. After fixing the pa- 
rameters (w_e,^ej0e) there still remain two impact pa- 
rameters. For e = and the chosen impact parameters 
iieiCe) the corresponding null geodesic arrives to J^^ at 
some moment of the retarded star time w+o in the di- 
rection {9o,4>o) with the impact parameters (^o,Co)- We 
focus our attention on the bending angles of a photon 
which are determined by (9o, 4>o)- For circularly polarized 
photons (e = ±|£| ^ 0) with the same initial scattering 
data the arrival angles (0(±), </'(±)) are slightly different, 

9(±)^9o±\e\M, + 0{e^), 

= 4>o± \e\^(t>o + O(e') . (94) 

The quantities A^o and A0o allow us to measure polar- 
ization shift of the bending angles. We shall analyze how 
the polarization shift of the bending angles depends on 
the impact parameters ($e,Ce)- 

Since for s 7^ we do not know analytical solutions, 
to solve this problem we use numerical methods. There 
are certain restrictions which have to be taken into ac- 
count in the numerical calculations. First of all, we must 
chose the initial impact parameters so that photons are 
not captured by the black hole. Note, that in the calcu- 
lation of the polarization contribution to the equation of 
motion the ultrastationary form of the metric was used 
[l|. For this reason, trajectory of a photon should not 
enter the ergosphere. Moreover, in the very vicinity of 
the ergosphere the effective force is growing fast, and 
the validity of the adapted linear in e approximation can 
be violated. To control the validity of the approximation 
we always check that the expansion ([M)) is valid and the 
deflection of the numerical result from the linear law is 
small. 

Before we describe the adopted numerical scheme, let 
us make a couple of additional remarks. We use the 



B. Approach and results 

The details of the scheme of the numerical calcula- 
tions are given in Appendix A. Here we briefly dis- 
cuss the main important steps of these calculations, the 
adopted choice of the initial data, and the obtained re- 
sults. We start the numerical integration by choos- 
ing the initial data for a null ray, that is, by fixing 
Ve = {w^eTUe = COS 9e , (jje , Pe, Xe} of au incoming at 
J^^ photon. Using the axial symmetry of the problem 
we put (f)e — 0- The Kerr metric is invariant under the 
reflection with respect to the equatorial plane y ^ —y. 
Equations (p5|) - ([7T|) also respect this symmetry, if one 
simultaneously changes e — ^ —e in these equations. For 
this reason, we consider only the values of 9^ lying in the 
upper semisphere < 9e < tt/2. To avoid the problem 
with bad behavior of the spherical coordinate (p near the 
symmetry axis, we always take a non-vanishing inclina- 
tion angle 9^- 

For the incoming null ray we use the advanced star 
time W- as a parameter along the ray. The inverse ra- 
dius X increases from its initial value a; = at to 
some value Xmax = M/rmin corresponding to the ra- 
dial turning point, where the radial coordinate reaches 
its minimal value rmin- After the turning point, we use 
the retarded star time w+ as a parameter along the ray 
(for details see Appendix A). In order to make w± well 
defined we chose the lower value in the integral in PSI) . 
which determines r^,, to be equal to 2M. Thus, at the 
turning point r = rmin we have 

w+ = {w_ - 2r,/M) . (97) 

At the end point at J^^ , where x again takes zero value, 
the retarded star time w+ determines the time of arrival 
of the null ray. In the coordinates ((59)) one has 

w+ 1 j+ = w- \j- + Aw- + Aw+ - 2 / -^dx . (98) 

Here Az«_ and Aw^ are calculated along the ingoing and 
outgoing null rays, respectively. 



10 




FIG. 3. Polarization shift of the bending angles. The left column contains plots showing dependence of the shift parameter 
Aa (shift angle divided by |e|) on the angle x for four different inclination angles 9^- The right column contains similar plots 
for the parameter /? G [0, 27r), which defines the direction of the shift. Figures (al) and (a2) correspond to 6e = tt/IO, Figs, 
(bl) and (b2) correspond to 6^ = n/6, Figs, (cl) and (c2) correspond to 9e — vr/3, and Figs, (dl) and (d2) correspond to 
9e = 7r/2. The vertical dashed lines in Figs, (dl) and (d2) cut out the values of x close to tt. 



The numerical procedure described in Appendix A al- 
lows one to calculate the polarization shift of the benging 
angles. Here we present only the results for the extremal 
Kerr spacetime (a = AI,a = 1), where the expected ef- 
fects are the most profound. To illustrate the dependence 
of the shift on the inclination angle 9e , we present the re- 
sults for its four different values, 9e = tt/IO, 7r/6, 7r/3, 
and 7r/2. 



For a chosen inclination angle 9,, one still needs to 
specify the impact parameters, which belong to the two- 
dimensional impact plane. In order to reduce the space 
of parameters we use the following approach. We fix 
the impact radius p and perform calculations for a dis- 
crete set of the angles x- By fixing p we keep the value 
of the asymptotic total angular momentum \L\ constant 
[see Eq. (pij) ]. Figure [3] shows the results obtained for a 
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FIG. 4. Capture domain for an extremal black hole (a — 1) 
and the inclination angle 9^ = n/2. The circle defines the 
'tube' of the null curves with the impact radius p = 7. Incom- 
ing null geodesies with impact parameters within the shad- 
owed region (capture domain) are captured by the black hole, 
while for the impact parameters outside this region scattered 
photons reach J"^ . 



special choice of p = 7. The argument in the plots pre- 
sented in the Figure is the angular parameter x, which 
parametrises the 'tube' of null curves with fixed p. To 
better illustrate our results, we connected the points gen- 
erated by the numerical computation into the smooth pe- 
riodic in X € [0, 27r) curves constructed by the spline in- 
terpolation of the ninth degree and followed by the poly- 
nomial interpolation using the Lagrange method. The 
function /3 is a continuous function of x which takes its 
values on a circle 5*^. To avoid its discontinuity in the 
special case shown in Fig. (d2) for the inclination angle 
Be = 7r/2, we shifted its range from [0, 2tt) to [tt/2, 5tt/2). 

For a chosen inclination angle 6^ a photon trajectory 
is defined by the impact parameters (^e,Ce)- For small 
impact parameters the photon is captured by the black 
hole. This means that one can identify two regions on 
the impact plane with qualitatively different photon tra- 
jectories. We call the first region, for which the capture 
occurs, a capture domain. It is easy to show that by the 
reflection t ~t,(p ^ —cj) this region is connected with a 
region on the outgoing impact plane, which is called black 
hole shadow (see, e.g., [23], [1^, and [IJl). The capture 
domain for the inclination angle 9^ = 7r/2 is shown in 
Figure S) This is the case when the incoming photons 
begin their motion parallel to the equatorial plane with 
(j)e = 0. The circle in this picture defines the 'tube' of the 
null curves with the impact radius p = 7. It touches the 
capture domain at the point where x = ^- Null geodesies 
emitted with the impact parameters close to this point 
revolve around the black hole many times before they ei- 
ther get captured or reach a remote observer. This means 
that for the chosen values of p = 7 and x close to tt the 
bending angles become large. In plots (dl) and (d2) we 



cut out the region of x near tt, where the bending angles 
become larger than tt. One can see that the correspond- 
ing polarization shift Aa grows fast near x — n (see plot 
(dl) in Figure [31). However, our calculations show that 
inside this region Aa can reach the value of ~ 5.6 still 
preserving the linearity in e with the accuracy w 97.5%. 

It is possible to show that for other values of the incli- 
nation angle 9e the impact radius p corresponding to the 
points inside the corresponding capture regions is always 
smaller than 7. For the data presented in Figure [3] we 
kept the same value p = 7 and changed the inclination an- 
gle. All the plots are periodic in the impact angle x with 
the period 27r. When the inclination angle changes from 
7r/2 to the plots for Aa become more smooth and the 
maximal value of Aa decreases. For x — '^/'^ and 37r/2 
the azimuthal angular momentum vanishes, i.e., p — 0. 
Null geodesies with such value of p pass through the axis 
of symmetry [cf. Eqs. ([75)) and ([7^ ]. To avoid the prob- 
lem of the coordinate singularity of (/), we performed the 
calculations only in the vicinity of these points. 

Plots shown on the right side of Figure [3] illustrate the 
dependence of the angle /3 on the impact angle x- One 
can see that the range of /3 grows when one considers 
inclination angles close to 7r/2 (equatorial plane). 

As it was mentioned in the Introduction, Mashhoon 
estimated the polarization shift of the bending angles in 
the weak field approximation [3-5]. In our notations, the 
estimated relation (see Eq. ([Ij) can be written in the 
following form: 

1 A I I 47ra|£| 

separation angle ~ 2Aa £ ^ k ^ — , (99) 

p-^ 

where k is a dimensionless parameter. For an extremely 
rotating black hole (a = 1) the parameter k can be de- 
fined as 

.^^^ (100, 

For the chosen value p = 7 one has k « 55Aa. Ac- 
cording to the data presented in Figure [31 Aa changes 
approximately in the range 0.01 — 0.06, so that the cor- 
responding value of K is roughly between 0.6 and 3.3. 
The increase of the parameter k can be explained by the 
fact that the rays with larger value of Aa come closer to 
the black hole. 



VI. POLARIZATION SPLITTING OF IMAGES 
AND ARRIVAL TIME 

A. Image of a point-like source 

Let us now discuss a different problem. Let us again 
consider a spacetime with a rotating black hole. Suppose 
now that there is a bright point-like object emitting light 
isotopically and an observer which registers that light. 
For example, the observer takes a picture of this object 
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by using a telescope. If the light passes close to the black 
hole, position of the image will depend on impact pa- 
rameters of the corresponding ray. Suppose now that the 
emitted light is circularly polarized. Since trajectories 
of circularly polarized photons are slightly different than 
those of photons following null geodesies, they arrive to 
the observer at slightly different angles. As a result, the 
position of the corresponding images in the picture will 
be slightly different. 

Let us assume that the distances from the black hole 
to the emitter and to the observer are large. In this case, 
in order to describe the polarization dependent splitting 
of the images it is convenient to use the impact param- 
eters (^o, Co) which determine the position of the image. 
In this idealized set-up the emitter is located at . As 
before, we put (/)e = 0, so that its position is character- 
ized only by the inclination angle Q^. Denote the angles 
which determine direction to the observer by (0o, '/'o)- In 
general, there exist many null geodesies connecting the 
emitter and the observer. They differ by a winding num- 
ber describing how many turns around the black hole a 
photon makes before reaching the observer. Brightness of 
the photon images decreases with increase of the winding 
number (see, e.g., HJ]). We shall focus on one of these 
null geodesies. In this case, the initial and the final im- 
pact parameters are uniquely defined. These parameters 
determine the asymptotic integrals of motion (pe, q^) and 
(Po, 9o) [see Eqs. (PTj) and ([Ml) ]. For a null geodesic of the 
Kerr spacetime these quantities are integrals of motion, 
so that po = and qo = Qe- 

The situation is different for motion of circularly polar- 
ized photons. Because of the axial symmetry, p remains 
an integral of motion, whileg is not, qo qe- To analyze 
how the parameter q changes, we introduce the following 
function q{£): 

2 2 

, P y 2 2 /I ni N 

+ Tj - a y ■ (101 

(1 - r) 

This function is defined by solving the relations (|76p and 
()77|) for q. Certainly, to obtain a similar expression for 
q{£) one can use the other pair of equations, Eqs. ([75|) 
and ([77)). However, the constraint (I55|) guarantees that 
these expressions for q{£) are identical. The function 
(jlOip coincides with q^ and qo at J'~ and , respec- 
tively. 



B. Calculation scheme and results 

The calculation of the polarization dependent splitting 
of images is a two-point boundary value problem. In our 
calculations we proceed as follows. First we chose a null 
geodesic connecting the points {Oe,(t>e = 0) and {Oo,(t>o) 
and determine the corresponding values of the parame- 
ters {pe,qe), which are constant along the geodesic. To 



find a null ray connecting the same points for e 7^ 
we use the shooting method and implement a numerical 
code written in Fortran (see, e.g., [2^). Namely, we use 
the initial seed values of Pe and qe and using the proce- 
dure described in Appendix A integrate the dynamical 
equation ()69|) up to the final point. We compare the nu- 
merical result for the final angles (9, (f>) with the required 
final data (9o,(t>o)- This comparison allows us to correct 
the seed values of Pe and qe and repeat the integration, 
until the required accuracy is achieved (for more details 
of the shooting method see, e.g., Chapter 17 of the book 
[29! ) . The numerical integration of the dynamical equa- 
tion uses the fifth-order Runge-Kutta method with the 
accuracy of 10^^. 

As earlier, we consider an extremal black hole, a = 1. 
To illustrate the results of the calculations, we consider 
a special case when the inclination angle for incoming 
photons is 7r/4, so that (j/e — cos9e, (j)e) = (a/2/2, 0). We 
also take e = ±0.01. The polarization splitting of the im- 
ages for different relative positions of the emitter and the 
observer is illustrated in Figure O The plots show posi- 
tion of the images in the impact plane for the following 
positions of the observer (j/o = cos9o,4'o)' (0.85,7.0), 
(0.4,4.5), (0.05,4.0), and (-0.12,3.8). The calculated 
values of the parameters (p(,,qe) are given in Table HI 
The splitting of the images can be characterized by the 
displacement parameter 

Ao(±) = |£rV(A^o(±))2 + (AC.(±))2, (102) 

where 

ACo(±) = Co(±) - "^0(0) 7 ACo(±) = Co(±) - Co(0) • (103) 

Here £,o(±) and Co(±) are the impact parameters of the 
left (— ) and right (-|-) circularly polarized photons, fo(o) 
and Co(o) correspond to the null geodesic. The values of 
the displacement parameter can be found in Table El 

Let us discuss photon trajectories for the same posi- 
tion of the emitter {\/2/2, 0) and for the case when the 
observer is located at the point (0.85, 7.0). We chose this 
case only for an illustration. Let us emphasize that in 
this case, the bending angle is larger than tt, that is, the 
photon passes the black hole at a relatively short dis- 
tance. This makes the polarization splitting of images 
more pronounced. However, the linear in e approxima- 
tion remains valid. Plots presented in Figures |BH7] show 
the photon coordinates x and y as functions of the an- 
gle (j). In Figure m the lower solid line demonstrates how 
photons starting at {x = 0) come close to the black 
hole, reach their turning point, where x takes its maxi- 
mal value, and propagate to J7^+ (x — 0). The upper solid 
curve shows position of the ergosphere with respect to a 
photon. To define the position of the ergosphere one uses 
Eq. where y — y{(j)) is defined by the photon tra- 

jectory. We calculated photon trajectories for a geodesic 
case (curve 0), and for the left (e = —0.01) and the right 
(e = -1-0.01) polarized photons, (curves —1 and -1-1). The 
difference between these curves is small. To demonstrate 
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-4.0566- 
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-2.8521 



-2.8517 
(11) 



-2.8513 



-6.2264 
-6.2265 
-6.2266- 
-6.2267 
-6.2268 
-6.2269 



+1; 



-3.5076 



-3.5074 
(IV) 



-3.5072 



FIG. 5. Images of a point-like emitter. Points —1 and -1-1 correspond to the left and right polarized photons, and point 
corresponds to a null geodesic. Figures correspond to the same initial direction of the emitter, {ye — cosde,cj)e) ~ (\/2/2, 0) 
and different positions of the observer {yo = cos6o,4>o)'- (0.85,7.0), ayo = —1 for Fig. (I); (0.4,4.5), ayo — 1 for Fig. (II); 
(0.05,4.0), aye = 1 for Fig. (Ill); (-0.12,3.8), ayo = 1 for Fig. (IV). 



TABLE I. Numerical data for Fig. El 



Point 


e 


Pe, plot (I) 


ge, plot (I) 


Pe, plot (II) 


(je, plot (II) 


Pe, plot (III) 


ge, plot (III) 


Pe, plot (IV) 


ge, plot (IV) 


-1 


+ 1 


-0.01 


-F-0.01 


1.978342 
1.980668 
1.982997 


11.739188 
11.721593 
11.703334 


2.614062 
2.613719 
2.613382 


17.598726 
17.593052 
17.587374 


3.086741 
3.086517 
3.086292 


28.016331 
28.012854 
28.009333 


3.482207 
3.482061 
3.481913 


38.936516 
38.933788 
38.931049 



TABLE II. Values of the displacement parameter Ao(±). 



Plot 






(I) 


1.4607 


1.4831 


(11) 


0.0756 


0.0754 


(III) 


0.0398 


0.0402 


(IV) 


0.0263 


0.0265 



formation about y = y{<t')- The numerical data for the 
asymptotic parameters (pe i Qe ) £^re given in Table IIIII 

TABLE III. Numerical data for Figs. [MSl 

Curve e pe ge 

-1 -0.01 1.978342 11.739188 
1.980668 11.721593 
+1 -1-0.01 1.982997 11.703334 



this, we presented a zoomed part of the trajectories near As we explained earlier, p is an exact integral of mo- 
their radial turning point. Figure [7] contains similar in- tion, while q changes along trajectory of a circularly po- 
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FIG. 6. Plot X = x{(j)). Curves —1 and +1 illustrate tra- 
jectories of the left and right circularly polarized photons, 
respectively. Curve illustrates a null geodesic. 



FIG. 8. Plot q — Qe versus (p. Curves —1 and +1 correspond 
to the left and right circularly polarized photons, respectively. 
Curve corresponds to a null geodesic. 




1 2 3 4 5 6 7 



FIG. 7. Plot y = 2/(0). Curves —1 and +1 illustrate tra- 
jectories of the left and right circularly polarized photons, 
respectively. Curve illustrates a null geodesic. 



larized photon [see Eq. (|10ip ]. Figure H] illustrates this 
change. To estimate an error in the linear in e approxi- 
mation we introduce the following quantities: 

Ap(±) =p(±|e|)-p(0), A<z(±) =(z(±|£|)-g(0), (104) 

and define the relative errors of the first order expansion 
in e 

_ |Ap(+) + Ap(_)| _ |Ag(+) + Ag(„)| 

'^<±'- 2|Ap(±)| ' '^<±'- 2|A9(±)| • 

(105) 

Using the data presented in Table UlIl we derive Sp^_^ « 
0.07%, <5p,^, « 0.06%, (5,(_) « 1.89%, and (5,,^, « 1.82%. 
This means that non-linear effects are small and can be 
neglected. 

C. Arrival time 

Let us briefly discuss the last of the three problems 
formulated at the end of Sec. IV. Namely, how the ar- 



rival time depends on the photon polarization. In order 
to determine this quantity, one can use the relation (|98|) , 
where the integration is performed along a photon trajec- 
tory. If one parameterizes the photon trajectory by the 
star time parameters, as it is done in our calculations 
(see Appendix A), the quantity 

Aw ^ w-^-\J+ — w-\j- (106) 

is obtained simply as the difference between the initial 
and the final values of the corresponding star time pa- 
rameters. Let us again emphasize that the value of Aw 
depends on the choice of the upper limit in the integral in 
([M]) . which we chose to be equal to 1/2, that corresponds 
to r = 2M. However, when one compares the registration 
times of photons of different polarizations emitted simul- 
taneously, the difference in their arrival times is uniquely 
determined. One has 

At(±) = Aw(±) - Aw(o) . (107) 

Here we present the results of the time delay of the 
polarized photons for the one of the above cases, when 
the spin-optical effects are more pronounced. The related 
data are given in Table IIIIl As before, we use the sub- 
script (0) for a null geodesic, while the subscripts (— ) and 
(+) stand for the left and right polarization of photons, 
respectively. For the null geodesic, we have the following 
value of the parameter (|106p : 

Aw(o) = 13.458554. (108) 

For the above data one has 

Ar(_)/e2 sa 0.17, Ar(+)/e2 « 0.18 . (109) 

These relations imply that polarized photons arrive later 
than 'non-polarized' ones, and that the time delay is the 
effect of the second order in the parameter e. Both the 
conclusions are in agreement with the general results on 
the time delay discussed at the end of Sec. II. 
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VII. DISCUSSION 

In this paper we applied the modified geometric op- 
tics approach developed in [ij to the problem of scat- 
tering of polarized light by a rotating black hole. We 
demonstrated that trajectories of polarized photons are 
null curves, which coincide with null geodesies only in 
the limit a; — >■ cx). A deviation of the null rays from 
null geodesies is controlled by the small parameter e = 
±(2cL'Af)~^. Since the polarization dependent term in 
the photon propagation equations (|69I) - ([7T|) rapidly falls 
down at infinity, the initial I?e and the final T^o scat- 
tering data for null rays are the same as those for null 
geodesies. These data include the (null) time and direc- 
tion at , as well as a vector in the 2D space of impact 
parameters. Working in (3-1-1) formalism we studied the 
equations for trajectories of polarized photons in the Kerr 
geometry. We developed a solver, which allows one to in- 
tegrate the required equations and to solve the scattering 
problem, that is, to find numerically the scattering oper- 
ator ^'(e) : Ve ^ Vo- We used this approach to study 
how bending angles of the scattered photons depend on 
their polarization. We also analyzed how the splitting of 
the image of a point-like source depends on polarization 
and calculated the time delay for polarized photons. It 
should be emphasized that the photon polarization en- 
ters the dynamical equations through the factor e, which 
depends on the photon frequency. This means that all 
the above discussed eflfects are frequency dependent. 

Unfortunately, in the framework of the approach used 
in the paper, one cannot consider effects for polarized 
photons passing quite close to a rotating black hole, e.g., 
propagating through its ergosphere. This happens be- 
cause the calculations are based on the (3-f 1) split of the 
spacetime obtained by projecting it along Killing trajec- 
tories. This split does not work on the ergosurface and 
inside of it. In principle, one can avoid this problem, for 
example, by using another Killing vector ^ = ^(•(■j-t-rj^^^-j, 
with properly chosen constant 17, so that the new Killing 
vector is well defined inside the ergosphere. After that, 
one can glue the trajectories corresponding to the mod- 
ified GO approach obtained with the help of ^ and ^(j) 
in the spacetime domain where both of the vectors are 
timelike. Another, more interesting problem is to develop 
the modified GO approach directly in a four dimensional 
spacetime. Such approach would allow one to study the 
polarization dependent effects for photons which closely 
approach a rotating black hole. In particular, it is in- 
teresting to study how polarization of light modifies the 
black hole shadow. One can expect that the black hole 
shadow will in fact be smeared for the polarized light, 
and if the 'bright radiation' behind the rotating black 
hole is not monochromatic, the position of the shadow 
rim will depend on the frequency of the radiation. In 
other words, one would observe a rainbow effect for the 
black hole shadow. 

The polarization splitting of images might, in principle, 
be detectable by astronomical observations. However, 



this depends on sensitivity of the measuring devices. It is 
interesting to study possible applications of the obtained 
results to real astrophysical problems. 



Appendix A: Details of the numerical computation 

Here we present the main steps of the numerical com- 
putation. Equations (IM1) - ([7T|) determine the coordinates 
of a null curve {x, y, (j), t) as functions of the proper length 
£. Direct integration of the equations in this form is not 
convenient, since for the scattering problem £ changes in 
an infinite interval. To deal with finite quantities param- 
eterizing the entire null curve (from to J^~^), we use 
the star time coordinates w±. Namely, we use w_ for 
the part of the null curve from to its radial turning 
point, and w+ from this point to The condition (^7)) 
uniquely relates these parameters at the turning point. 

The complete set of the equations for x, y, and (f) as 
functions of w± written in the first order form is 

dx X dy Y 



dw± A ' dw± A 
dX 1 



dw± 2A-fj;, 



2-(p-7^)^-^^^(curl7), 



2A7..7^/- A^ 
^yy{cw[\ curl7) 



-Y. 



1 



dY 



dw± '^A'^y 



'{lyy,y^ '^^yy ,xXY ^xx,yX ) 
iP - 1<P) , 



^^^'^ (p-70)^ + ^^^r^(curl7). 



1xx{cmt:\ curl7)0 



-A, 



dw± Ajt/xj) 



where 



A=l 



70(p-7<^) V 



700 



x^6 



\X\. 



The constraint ((65l) . which takes the form 



-fx^X^ + JyyY^ + 



{P - l<t>? 



= 1, 



(Al) 



(A2) 



(A3) 



is used to control the accuracy of the calculations. 

In principle, one can use these equations to start the 
integration from a point with x = corresponding to 
J~ . In fact, we found that it is more convenient to use a 
slightly modified approach, which allows one to get higher 
precision. Namely, in the vicinity of Sf^ we parameterize 
the coordinates x and y by the angle cf), where for a ro- 
tating black hole this parameter is monotonic along null 
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curves. For this parametrization the initial values for a 
null curve at J~ are 



i=0 



dx 
dy 



0=0 



P 

p 



(A4) 



y/{q + a'^y'^){l - y'}) - p'^yl 



The asymptotic integrals of motion p and g, which enter 
these relations are related to the impact parameters by 
means of the expressions ((9T|) and (|93p . 
The equation for is of the form 



where 



A 



d(f) dcf) \ d(f) 



(A5) 



B = 

C = 

£ = 



700,: 



700 

27a; a; 
700,2: 



7:r 



2jxx (P ~ 70)\/7 

_ lxx.,y I lyyl4>4> 



(P- 70)^/7 



(curl7)j, , 
(curl 7):,; , 



(A6) 



7yy700 



(curl curl 7)^ 



272:2; (p-70)v^ 



(curl7)y. 



A similar equation can be written for 2/(0). This equation 
is obtained if one makes an interchange x ^ y in all the 



expressions that enter (jA5l) and (jA6l) and reverse the sign 
of the terms containing curl 7 and curl curl 7. To solve 
these equations we use the seventh-eighth order Runge- 
Kutta method with the accuracy of 10^^ and check our 
computation by using the constraint (IA3p adapted to this 
form of equations. We integrate the equations from the 
point where a: = up to the point where x = xi is of 
the order of lO^'^. Since one cannot guarantee that (j) is 
a monotonic coordinate along the entire null curve, we 
use the result obtained at the point Xi to find the initial 
conditions for the system (lAip . Then we integrate this 
system by using the seventh-eighth order Runge-Kutta 
method with the accuracy of 10~^ and check our com- 
putation by using the constraint (|A3|) . We integrate the 
equations through a radial turning point, until we reach 
the point where a; = a;2 is of the order of 10^^. At this 
point we derive new initial conditions and use them to 
continue the calculations by integrating the system (jASP 
up to the end point where x w 10^^. This allows us to ob- 
tain with the required accuracy the angular coordinates 
(6*0, <j)o) at ^7+ for the outgoing null rays. 
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